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We present experimental results attempting to fingerprint nonanalyticities in the magnetization 
curves of spin glasses found by Katzgraber et al. [Phys. Rev. Lett. 89, 257202 (2002)] via zero- 
temperature Monte Carlo simulations of the Edwards- Anderson Ising spin glass. Our results show 
that the singularities at zero temperature due to the reversal-field memory effect are washed out by 
the finite temperatures of the experiments. The data are analyzed via the first order reversal curve 
(FORC) magnetic fingerprinting method. The experimental results are supported by Monte Carlo 
simulations of the Edwards- Anderson Ising spin glass at finite temperatures which agree qualitatively 
very well with the experimental results. This suggests that the hysteretic behavior of real Ising spin- 
glass materials is well described by the Edwards- Anderson Ising spin glass. Furthermore, reversal- 
field memory is a purely zero-temperature effect. 



PACS numbers: 75.50.Lk, 75.40.Mg, 05.50.+q 



I. INTRODUCTION 

While the nonequilibrium behavior of spin glasses-!* 2 -' 3 - 
has been studied in detail and several aspects have been 
applied to fields as wide as biology and financial anal- 
ysis, less work^&ii 8 - has been done to understand the 
hysteretic behavior 9 of spin glasses in a time-dependent 
field. Because the understanding of hysteretic systems 
plays a prominent role in the development of magnetic 
recording media, understanding fundamental properties 
of simple hysteretic systems such as spin glasses is of 
paramount importance. 

In Ref. 8 Katzgraber et al. reported on a novel memory 
effect — the reversal-field memory effect — found in spin- 
glass systems: When the field is decreased from satura- 
tion to a reversal field Hr, upon return to saturation the 
magnetization exhibits a singularity (kink) at H = — Hr. 
The underlying spin-reversal symmetry of the Hamilto- 
nian is the source of this effect which was observed first 
numerically at zero temperature for the two-dimensional 
Edwards- Anderson Ising spin glass. 10 In contrast, sys- 
tems which do not possess spin-reversal symmetry, such 
as the random- field Ising model ^ do not exhibit this 
memory effect. Because the kinks in the magnetization 
curves are easily overlooked, Katzgraber et al. used the 
first order reversal curve (FORC) method 8 - ' 11 ' 12 to char- 
acterize the reversal- field memory effect. In a FORC dia- 
gram, the kinks in the magnetization curves are captured 
by a pronounced vertical ridge along the coercivity axis 
(see below). 

Using the superb fingerprinting abilities of the FORC 
method we have attempted to experimentally detect 
the reversal-field memory effect in Feo.sMno.sTiOa Ising 
spin-glass samples. Our finite-temperature experimen- 
tal results show that the reversal-field memory effect 
is not present. Moreover, these results are verified via 
finite-temperature nonequilibrium Monte Carlo simula- 



tions of the three-dimensional Edwards-Anderson Ising 
spin glass. The numerical results show that the reversal- 
field memory effect can only be observed at experimen- 
tally inaccessible temperatures. The good agreement be- 
tween the simulation and experiments suggest that the 
hysteretic properties of Feo.sMno.sTiOa can be well mod- 
eled using the short-range Edwards- Anderson Ising spin 
glass model. 

The paper is structured as follows. First we review the 
FORC method in Sec. [TTJ In Sec. IIIII we present experi- 
mental results on Feo.sMno.sTiOs. In Sec. IIVI results of 
zero- and finite-temperature Monte Carlo simulations of 
the three-dimensional Edwards- Anderson Ising spin glass 
are presented, followed by concluding remarks in Sec. fVl 



II. OUTLINE OF THE FORC METHOD 

FORC diagrams ; 8 ' 11 : 12 which can be viewed as 
a model-independent generalization of Preisach 
diagram a 13 ' 14 complement current methods to charac- 
terize magnetic interactions in hysteretic systems ) 15 ' 16 : 17 
such as the SM method J&l 9 - The advantage of the FORC 
method over other approaches is the extreme sensitivity 
to microscopic details of hysteretic systems. Although 
the wealth of information a FORC diagram delivers 
still remains to be fully characterized, the method can 
be used as an extremely sensitive "fingerprint" of a 
hysteretic system. 

To calculate a FORC diagram, a family of FORCs 
with different reversal fields is measured cither ex- 
perimentally or numerically^ The measurements start 
at saturation and the field is then reversed at different 
values of the reversal field Hr. The mixed second-order 
derivativ e) 11 ' 12 of the magnetization M(H, Hr) as a func- 
tion of the applied and reversal field yields the FORC 
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FIG. 1: (Color online) Zero field cooled (ZFC), field cooled 
(FC) and thermoremanent magnetization (TRM) of the 
Feo.sMno.sTiOa sample as a function of temperature T for 
an applied field of H = 50e. For this experimental sample 
the glass transition occurs at T c ~ 21 K. 

diagram p(H,H R ) given by 

p(H, H R ) = — i [d 2 M / dHdH R ] . (1) 

In general, a rotation of the coordinates to 

H c = [H-H R }/2, (2) 
H h = [H + H R }/2, (3) 

the local coercivity H c and bias i?b, respectively, yields 
the FORC distribution p(Hh,H c ). Therefore, under the 
assumption that a hysteretic system can be described 
by microscopic switching units, each of which have a 
given coercivity H c and a bias field ift>, the FORC di- 
agram corresponds to the distribution of these micro- 
scopic bias/coercivity fields. In the case of discrete 
spin models such as the random-field Ising modelj££ 
these microscopic switching units can, in general, be 
identified with the simulated spinsj^I This is not the 
case for systems which exhibit frustration, such as the 
Edwards- Anderson Ising spin-glass model. FORC dia- 
grams have been applied to a variety of hysteretic sys- 
tems and fields ranging from geological applications^ to 
magnetic materials ^ 22 ' 23 as well as more exotic systems 
such as cyclic voltammetry^ superconducting perovskite 
materials^ and thermal hysteresis^ 



III. EXPERIMENTAL RESULTS 

The experiments have been performed on a single 
crystal of the Ising spin glass Feo.sMng.sTiOs in a 
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FIG. 2: (Color online) Hysteresis curve of Feo.sMno.sTiOs at 
6.5 K, which corresponds to a temperature of T w 0.31T C (see 

Fig. Q5. 



Quantum Design MPMS-XL 5T SQUID (superconduct- 
ing quantum interference device) magnetometer. Low- 
field magnetization vs temperature curves in zero field 
cooled (ZFC), field cooled (FC), and thermoremanent 
magnetization (TRM) protocols for Feo.sMno.sTiOs are 
shown in Fig. [T] The spin-glass transition tempera- 
ture of Feo.5Mn .5Ti03 is T c « 21K. These measure- 
ments are made with the applied field parallel to the 
c axis of the sample — the corresponding curves mea- 
sured perpendicular to the field show only a weakly 
temperature-dependent paramagnetic response, no irre- 
versibility between the ZFC and FC magnetization and 
zero remanencei2£ In addition, in Fig. [5] a magnetization 
(M) vs applied field (H) hysteresis curve measured at 6.5 
K, the temperature of the FORC experiment described 
below, is shown. An extensive study of the hysteretic 
behavior of Feo.sMno.sTiOa has been reported by Ito et 
al. in Ref . [H. 

The FORC measurements have been made at as low 
a temperature as possible to minimize influences on the 
FORC distribution from thermal relaxation effects. The 
temperature must, however, be high enough for the sam- 
ple to have a reversible magnetization at the highest mag- 
netic field available (5 T). The chosen working tempera- 
ture of 6.5 K (T/T c w 0.31) optimally fulfills these two 
requirements. 

The derived FORC distribution is shown in Fig. [3] 
The individual FORCs behind the distribution have been 
measured using field steps of 80 mT and a limiting field 
of ±5 T. The measurements are made at a constant field 
increase rate governed by the time it takes to change the 
field and record one data point (each data point must 
be recorded at constant field in the SQUID). To avoid 
some of the influence of relaxation of the magnetization 
after decreasing the field from the initial high value to 
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FIG. 3: (Color online) FORC diagram of Feo.5Mno. 5 Ti0 3 at 
T — 6.5 K (T/T c w 0.31). Note that the distribution is not 
normalized. 



Hr, the sample is kept for 1000 s at the reversal field 
before the actual recording of data during field increase 
starts. A preliminary report on these results is found 
in Ref. [2^. It can be noted from Fig. [3] that there is a 
wide distribution of both coercivity and bias fields. The 
distribution shows a smooth variation with a weak max- 
imum at a coercivity of about 1 T and a bias field close 
to zero; and it only slowly decays to zero from a broad 
ridge along the coercivity axis. In particular, the horizon- 
tal ridge reminiscent of the reversal-field memory effect 
along the H c axis has "melted" in comparison to the zero- 
temperature results from Monte Carlo simulations, and 
the distribution has broadened along the vertical Hy, axis. 
By comparing to Monte Carlo simulations at finite tem- 
peratures (see below) we therefore conclude that finite 
temperatures and therefore fluctuations of the magnetic 
moments disrupt the reversal-field memory effect consid- 
erably even at temperatures considerably lower than the 
transition temperature of the glass phase. 
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FIG. 4: (Color online) FORC diagram of the three- 
dimensional Edwards- Anderson Ising spin glass at zero tem- 
perature (see Ref. |8| for further details) . The data show a pro- 
nounced ridge at zero bias (Hb = 0) due to the reversal-field 
memory effect which is a manifestation of the spin-reversal 
symmetry of the underlying Hamiltonian. 
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IV. NUMERICAL RESULTS 

The experimental results presented in Sec. IIHI do not 
agree with the numerical results presented in Ref. 8. In 
particular, the narrow ridge along the horizontal axis 
that captures the reversal-field memory of the Edwards- 
Anderson Ising spin glass is washed out (see Fig. HJ. To 
better understand these findings, we simulate the three- 
dimensional Edwards- Anderson Ising spin glas o 1 ' 10 given 
by the Hamiltonian 

H=J2JioSiS j -Hj2Si (4) 

(id) » 

at finite temperature. Here Si £ {±1} represent Ising 
spins on a cubic lattice of size N — L 3 with periodic 



FIG. 5: (Color online) FORC diagram of the three- 
dimensional Edwards- Anderson Ising spin glass at T = 0.30, 
which corresponds to T/T c m 0.31. The sharp ridge found at 
zero temperature (Fig. 2} is washed out completely yet the 
asymmetry along the horizontal axis remains, in very good 
agreement with the experimental data shown in Fig. [3] 



boundary conditions. The nearest-neighbor interactions 
between the spins Jij are chosen from a Gaussian dis- 
order distribution with zero mean and standard devia- 
tion unity. H represents an externally applied field. The 
model has a spin-glass transition at T c « 0.95^ 

The zero-temperature dynamics of the Edwards- 
Anderson Ising spin glass is simulated by changing the 
external field H in small steps starting from positive sat- 
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uration. After each field step we compute the local field 
hi of each spin 

//, ^./,,.S' ; //. (5) 

3 

A spin is unstable if it points opposite to its local field, 
i.e., if hiSi < 0. Randomly chosen unstable spins are 
flipped and their local fields updated until all spins are 
stabled 

At finite temperatures the system is simulated with a 
Monte Carlo approach ! 32 ' 33 the external field is changed 
in small steps and for each field step the system is equili- 
brated using heat-bath Monte Carlo. Tests show that for 
1000 Monte Carlo lattice sweeps at T = 0.30 we obtain 
hysteresis loops which are independent of Monte Carlo 
time within error bars. While the system is not in full 
equilibrium 34 neither are the experiments and so we ar- 
gue that our results describe the intrinsic noncquilibrium 
nature of the finite-temperature experiments correctly. 
For both zero and finite temperatures we choose a sat- 
uration field f/gat = 16.0 and perform 400 field steps. 
The presented data for 20 3 spins are averaged over 2000 
disorder realizations. 

Figure 2] shows a numerical FORC diagram at zero 
temperature for the Edwards- Anderson Ising spin glass 
(see Ref . H) . The data show a pronounced ridge at = 
which is reminiscent of the reversal-field memory effect. 
In Fig. [5] we present data at T = 0.30, a temperature 
(in dimensionless units) which agrees with the values 
used in the experiments presented in Sec. IIIII The sharp 
horizontal ridge is completely washed out in very good 
qualitative agreement with the experimental results and 
shows that for the Edwards-Anderson Ising spin glass 
reversal-field memory is destroyed by finite-temperature 
fluctuations. This suggests that the reversal-field mem- 
ory effect 8 can only be probed at temperatures much 
lower than can currently be achieved experimentally. 

V. CONCLUSIONS 

The hysteresis behavior of spin glasses is governed by 
temperature-relaxation phenomena and field sweep rates 



dictate the width and shape of the hysteresis loop at fi- 
nite temperatures. We find fundamental differences be- 
tween the field-driven hysteresis derived at zero temper- 
ature by Monte Carlo simulations and the experimen- 
tally as well as numerically observed behavior at finite 
temperature. In particular, the FORC distributions of 
finite-temperature Monte Carlo simulations and experi- 
ments on Feo.sMno.sTiOa (an Ising spin glass material) 
at T/T c w 0.31 agree qualitatively well (see Figs. [5] and 
[3]). The zero-temperature behavior which is dominated 
by the reversal-field memory effect and characterized by 
a sharp ridge along the horizontal axis of the FORC di- 
agram is absent at finite temperatures. 

Our results therefore show that the reversal-field mem- 
ory effect can only be observed at temperatures close 
to or equal to zero which are inaccessible experimen- 
tally with current technology. Furthermore, the quali- 
tative agreement between the finite-temperature Monte 
Carlo data and the experimental results suggest that 
Feo.sMno.sTiOa is well described by a (short-range) 
three-dimensional Edwards- Anderson Ising spin glass. It 
would be of interest to further characterize other mate- 
rials as well as effective models attempting to describe 
these using the FORC method. This is of paramount im- 
portance for materials for which model Hamiltonians are 
currently either unknown or under debate. For example, 
Feo.85Zno.15F2 is a diluted antiferromagnet in a fiel d 35 ' 36 
which is expected to be well described by a random-field 
Ising modeli^' 6 - Currently experiments as well as simu- 
lations are being performed to characterize this mate- 
rial/model using the FORC method. 
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